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Abstract 

The total cross section for electron-positron annihilation into hadrons is calculated in the region 
just below the B meson threshold. QCD corrections up to third order, quark mass effects, initial 
state radiation, the running of chqed an d the tails of the T resonances are included in the prediction. 
For given a s the prediction is accurate to 0.5 percent. An experimental measurement with the 
corresponding precision would allow to determine a s with high accuracy at intermediate energies. 



1. Introduction 

One of the most precise and from the theoretical point unambiguous determinations of the strong 
coupling constant a s has been obtained from the hadronic decay rate of the Z boson. Ingredients are 
the large event rate which allows to reduce the statistical uncertainty in a s down to 0.003, the precise 
calibration through the accurate measurement of the luminosity or the leptonic rate jO and, last but 
not least, the elaborate theoretical calculations ||. This high energy determination of a s should 
and could be complemented by a similar measurement of -Rhad i n the energy region just below the B 
meson threshold: large event rates are available at the e + e~ storage ring CESR and an experimental 
analysis is currently underway Q. Theoretical predictions for i?had below the BB threshold have 
been presented earlier [||. These include terms up to order for the massless limit and for m 2 /s 
corrections as well, terms of order a 2 s for the quartic mass terms (m 2 /s) 2 , terms of order a s (m 2 /s) 3 and 
contributions from virtual bottom quarks up to order a 2 (s/m 2 ). Recently it has been demonstrated 
that this expansion in (m 2 /s) provides an excellent approximation not only around -^/s = 10 GeV, 
but also down to 6 or even 5 GeV. 



The forthcoming experimental analysis necessitates on the one hand a thorough understanding 
of the background, e.g. from two photon events or the feedthrough from r pairs into the hadronic 
cross section, and, on the other hand, a realistic prediction for the annihilation channel. The former 
is outside the scope of this paper, the latter is the subject of this work. It requires to incorporate 
the effect of the running QED coupling constant and initial state radiation. The photon vacuum 
polarization leads to an increase of the cross section by about 7-8%. Magnitude and even sign of 
initial state corrections depend on the experimental procedure, in particular on the minimal mass 
of the hadronic system accepted for the event sample. Another important issue is the treatment 
of the tails of the T resonances, which contribute coherently through mixing with the photon and 
incoherently through their radiative tails. In order to study these effects and the associated theoretical 
uncertainties separately the corresponding analysis for the reaction e + e~ — > t + t~ is presented in 
section 2. These results may also serve as an independent experimental calibration of the cross 
section. The corresponding predictions for the hadronic cross section are discussed in section 3. In 
section 4 we comment on the accuracy of the presented approach and section 5 contains our summary 
and conclusions. 



2. Lepton Pair Production 

To calibrate the predictions for the hadron production cross section it seems appropriate and useful to 
calculate in a first step the cross section for lepton pair production. To arrive at a reliable prediction, 
initial and final state radiation, the effect of vacuum polarization and the influence of the nearby T 
resonances must be included. In this section results are presented for r pair production. 

(a) Initial state radiation 

The most important correction to the total cross section is introduced by initial state radiation. It 
leads to a reduction of the invariant mass of the produced lepton pair or hadronic system and, for 
fairly loose cuts, to a significant enhancement of the cross section, albeit with events of significantly 
lower invariant mass of the system of interest. For the precise determination of i?had discussed below 
it is advisable to exclude the bulk of these low mass events. This reduces the size of the correction 
and at the same time the dependence on the input for i?had from the lower energy region. 

The treatment of initial state radiation has advanced significantly as a consequence of the 
detailed calculations performed for the analysis of the Z line shape. The result for the inclusive cross 
section can be written in the form 

o-(s) = f 1 dza (sz)G(z) . (1) 

Jzo 

The cross section including photon vacuum polarization (running a) is denoted by a$. The invariant 
mass of the produced fermion pair is given by s z, where 

^SSin <z <Z<l (2) 

s 

and m m i n = 2mg for lepton pairs and m m \ n = 2m w for hadron production. In the cases of interest for 
this paper zq will have to be adopted to the experimental setup. Typically it is significantly larger than 
the theoretically allowed minimal value. For E cm = 10.52 GeV a cut in z around 0.25 corresponding 
to roughly 5 GeV in the minimal mass will in the case of hadronic final states exclude the charmonium 
resonance region and the charm threshold region as well, thus limiting the hadron analysis to truely 
multihadronic final states. 



The complete radiator function G(z) up to order a 2 has been calculated in || [7j. The resum- 
mation of leading logarithms is discussed in || (see also ||] ) . For the present purpose an approximation 
is adequate which is exact in order a and which includes the dominant terms of order a 2 plus leading 



logarithms. For the radiator function G we thus take [1C 
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where 



r(i + /3)' 

0.5772 ... is Euler's constant. Eq. 



is suited for quick numerical integration. The 
difference to the initial state radiation convolution using the complete order a 2 result (eq (3.12) of 
§§) is below one permille. The predictions for the reference energy of 10.52 GeV and a variety of cuts 
zq are listed in Table |]. The differential distribution da/dz and the cutoff dependent cross section 



Table 1: Dependence of the cross section e + e — * t + t on the cutoff zq for E\ = 10.52 GeV. 
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0.7831 


0.7359 
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0.9988 
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a(z > zq) are displayed in Fig. |T|a, b, respectively, as functions of z, zq. For the cross section before 
convolution the corrections from the vacuum polarization, the tails of the T resonances and final 
state radiation are included. Initial state radiation of lepton pairs and of hadrons has been treated 
in 0, |Il]| . These could easily be incorporated in the present formalism. The resulting corrections are 
at the permille level and will be ignored for the moment. Corrections from initial state radiation are 
evidently huge for loose or extremely tight cutoffs, with opposite sign. The strong cutoff dependence of 
the result is apparent from Table [l] and Fig. [l] as well. The result is relatively stable against variations 
of zq for values of zq below 0.75. At the same time the sensitivity to the input for <jq from smaller 
s' = s z is eliminated for this choice. From this viewpoint a value around 0.7 to 0.8 seems optimal for 
the experimental determination of i?had- On the other hand, if one allows to use input from smaller 
s' from other experiments or from the measured events with photons from initial state radiation, then 
also lower values of zq down to about 0.25 are equally acceptable. 

The distribution da/dz is flat in the region 0.25 < z < 0.75 and numerically small. An 
experimental analysis based on a cutoff zq = s'/s with a reasonably symmetric resolution in s' is thus 
insensitive towards the details of the resolution function. In fact, even an uncontrolled shift of the 
central value by Azq = 0.05 would lead to a 1% deviation of the cross section only. 
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Figure 1: (a) The differential cross section d<r(e + e - — * r + T~)/dz and (b) the cutoff dependent cross 
section o~(z > Zq) in nb for the centre of mass energy E\ = 10.52 GeV as functions of z and zq, 
respectively. 



(b) Final state radiation 

Initial state radiation is strongly enhanced, a consequence of the inherently large logarithm in the 
correction function (cf. eq. @). Final state radiation with loose cuts, in contrast, is typically of order 
a/ tt without a large logarithm. In the massless limit and without a cutoff on the invariant fermion 
pair mass the correction factor (1 + 3/4a/7r) amounts to two permille only. It depends, however, on 
the lepton mass and the cutoff. For the relatively loose cuts with zq between 0.25 and 0.8 one may 
well use the totally inclusive correction function n 1 ' defined by 
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In the case of hadron production, the main motivation of this investigation, photons from final state 
radiation will in general anyhow be included in the hadronic invariant mass. 



(c) Leptonic and hadronic vacuum polarization 

The leptonic vacuum polarization in one loop approximation is given by 

n 77 ( s ) := n 77 ( s ) - n 77 (o) = Y. N cQ)P{s,™f) , (8) 

/ 

with 

1 2m}\ ( n oi J3-1 



P(», m/) = _ n + ^M 2 + /31n £_). (9 ) 



For \s\ 3> mj the function P(s,m/) is well approximated by 



f , as y mp.(s,m / ) = -^ + ln^--^ + ie^. (10) 

The hadronic vacuum polarization is obtained from it's imaginary part via a dispersion relation 

~ 2 aq 2 f°° iWgO , / ni s 

To avoid complications that arise from the numerical evaluation of the dispersion integral over the 
data in the timelike region, the integral is evaluated at the corresponding value of s in the spacelike 
region. We use parametrisations from the evaluation of eq. (jll|) in the spacelike region as provided 



by 1 12, 13]. The results for the complete cross sections calculated with the different parametrisations 
from [O] and |l^] agree to better than 10 -4 in the region of interest. The error from the identification 
of spacelike and timelike ii^ad (q 2 ) should be small except for the 6 quark contribution, where the 
threshold is fairly close to the q 2 values of interest. Therefore we subtract the perturbative b quark 
contribution evaluated for spacelike q 2 and add the corresponding value for timelike q 2 = s. We have 
checked that this ansatz is in excellent agreement with the numerical evaluationQ of (|IT]) in the timelike 
region. Effects due to T resonances are described in detail below. 

The running a is then obtained from the real part of IT through 

a(s) = (12) 

1 - ReII(s) 

The relative shift in a(s) from the hadronic plus leptonic vacuum polarization is shown in Fig. [2] as a 
function of yfs. The size of the individual contributions is listed in Table § for the reference energy 
E\ = 10.52 GeV and for a few selected lower energies. The uncertainty in the cross section from our 
treatment of the hadronic vacuum polarization in the timelike region is estimated to be below two 
permille. 

(d) Narrow T resonances 

Up to this point we have not included the contributions from the T resonances. As a consequence 
of their close proximity their effect will be enhanced and thus should be discussed separately. The 
BreitWVigner amplitudes from the T resonances above and below the BB threshold interfere with the 
virtual photon amplitude and thus enhance or decrease the cross section by 

Kn 6Te M\s-M 2 ) 



* We thank H. Burkhardt for providing the numerical evaluation of ([11]) required for this comparison 



Table 2: Individual contributions to Rell(s) • 10 2 . 
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Figure 2: The relative shift in a{s) as a function of yfs from the hadronic plus leptonic vacuum 
polarization as described in the text. 



The individual contributions are listed in Table || for the two cms energies E\ = 10.52 GeV and 
E2 = 9.98 GeV after the convolution with initial state radiation. For E\ these interference terms 
happen to cancel to a large extent, leaving a minute correction with a relative magnitude around 
10 -4 . For E2, however, due to the close proximity of T(25) one receives a non-negligible negative 
contribution of about —0.4%. This effect could become relevant in precision tests. 

Also the radiative tails of the resonances with M < yfs have to be taken into consideration. 
These can be easily calculated with the radiator function eq. (||), taking either the full Breit-Wigner 
resonance 
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or the narrow width approximation 
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Table 3: Relative contributions from the T resonances to the leptonic cross section a(e + e — ► r + r ) 
for the energies E\ — 10.52 GeV and E 2 — 9.98 GeV and two different values of the cutoff T7^ m i n . 
Interference terms and radiative tails are listed separately. Also given are the continuum contributions 
and the resulting predictions for the total cross section. 





T(l) 


T(2) 


T(3) 


T(4) 


T(5) 


T(6) 


M [GeV] 


9.460 


10.023 


10.355 


10.580 


10.865 


11.019 


r e [keV] 


1.32 


0.576 


0.476 


0.24 


0.31 


0.13 


r to t [MeV] 


0.0525 


0.044 


0.0263 


23.8 


110 


79 


E 1 = 10.52 GeV, m min = 2m T = 3.554 GeV: 


Continuum: 


0.9256 nb 


Interf.: 


2.7- HT 4 


2.9 • 10" 4 


7.5 • 10~ 4 


-1.04- 10" 3 


-3.0 • 10- 4 


-9.7- 10~ 5 


J2 interf.: 


-1.3- 10- 4 


Rad. tails : 


5.5 • HT 4 


2.5 • 10" 4 


8.1 • 10" 4 








ra d. tails: 


1.61 • 10" 3 


Total: 


0.9269 nb 


E 1 = 10.52 GeV, m min = 5 GeV: 


Continuum: 


0.9033 nb 


Interf.: 


2.9 • HT 4 


3.1 • 10~ 4 


7.7 • 10- 4 


-1.07- 10" 3 


-3.1 • 10" 4 


-9.8- 10" 5 


J2 interf.: 


-9.7- 10" 5 


Rad. tails : 


5.7- 10~ 4 


2.6 • 10" 4 


8.3 • 10" 4 








rad. tails: 


1.65 • 10" 3 


Total: 


0.9047 nb 


E 2 = 9.98 GeV, m min = 2m T = 3.554 GeV: 


Continuum: 


1.0200 nb 


Interf.: 


6.4- 10~ 4 


-3.51 • 10" 3 


-4.5 • 10" 4 


-1.6 • 10" 4 


-1.6 • 10" 4 


-5.9- 10" 5 


J2 interf.: 


-3.70 • 10" 3 


Rad. tails: 


1.08 • 10~ 3 












rad. tails: 


1.08 • 10" 3 


Total: 


1.0173 nb 


E 2 = 9.98 GeV, m min = 5 GeV: 


Continuum: 


0.9947 nb 


Interf.: 


6.8 • 10" 4 


-3.60 • 10" 3 


-4.6 • 10" 4 


-1.6 • 10" 4 


-1.6 • 10" 4 


-5.8- 10" 5 


J2 interf.: 


-3.75 • 10- 3 


Rad. tails: 


1.11 • 10~ 3 












rad. tails: 


1.11 • 10" 3 


Total: 


0.9920 nb 



In total they raise the lepton pair cross section by about 0.16% or 0.11% for E\ and E 2 , respectively. At 
this point a brief comment ought to be made concerning the treatment of T(45). The approximation 
of neglecting the energy dependence of the width of the resonance is justified even for Ei, since the 
nominal width T = 24 MeV is significantly smaller than the difference between the cms energy E\ and 
the mass of the resonance: AE = 60 MeV. 



3. The Total Cross Section for Hadron Production 



As stated already in the introduction the emphasis of this work is on the energy just below the 
bottom meson threshold. Virtual b quark loops can be easily taken into account, and in fact it has 
been demonstrated that the hard mass expansion works surprisingly well for these diagrams — not 
only far below but even close to the nominal b quark threshold. On the other hand this energy is 
sufficiently far above the open charm threshold such that charm quark mass effects can be included 
through an expansion in powers of m c , if quadratic and quartic terms are incorporated. With this 
method the production cross section can be predicted reliably not only in the high energy region 
but even relatively close to threshold through an expansion in m 2 /s. This approach was suggested 
originally in fllj, |!|. The calculation of i?had to second order, including the full mass dependence 
fj~5|l has demonstrated the nearly perfect agreement between approximate and exact result for the 
coefficient of the a 2 term for E cm above 4m. In fact, even for E cm around 3m, which is around 
the lowest advisable value of the cutoff m m \ n = 5 GeV, the deviation of the a 2 coefficient from the 
complete result leads to a difference of less than five percent in the rate. For the energy around 10 
GeV this region contributes through the radiative tail only, and the approximation is thus adequate 
throughout. In total we thus use the following individual contributions 
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The formulae are evaluated for n/ = 4 with a s and the charm quark mass interpreted accordingly. For 
the massless case and the m 2 / 's terms the results are available in third order, for the quartic and m 6 /s 3 
terms in second and first order a s , respectively. Tables which list the numerical values of the running 
quark masses and the magnitude of the individual contributions can be found in |5], |3| together with 
the details of the matching between the theories with np = 4 and 5. 

The predictions for i?had as a function of E cra are shown in Fig. ^ for three values of the strong 
coupling constant, a s (M^) = 0.115, 0.120 and 0.125. Fig. ||a displays the contributions induced by 
the light plus charm quark currents (as defined in eq. (|i~6|)), whereas in Fig. ||b only the light (u, d and 
s) quark current contributions are shown. Note, that this figure by definition does not contain QED- 
corrections from initial state radiation and vacuum polarization, but includes the tiny singlet terms, 
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Figure 3: (a) Rh&d( s ) as defined in eq. (jiqj, and (b) only the contributions from the light (u, d and s) 
quark currents for different values of a s . 



which cannot be attributed to one individual quark species and the final state QED-corrections <5qed- 
The quark masses are chosen to be rn c (jn c ) = 1.24 GeV (n/ = 4) and ffib(jnb) = 4.1 GeV (n/ = 5) 

(rif=5) 

corresponding to pole masses of 1.6 GeV and 4.7 GeV, respectively, if al 1 ' (M z ) = 0.120 is adopted. 
A variation of the charm quark mass around the default value by +300 MeV/— 300 MeV changes -Rhad 
only by +0.006/ - 0.004 for y/S = 10.52 GeV and by +0.031/ - 0.028 for ^/s = 5 GeV. 

Throughout this paper the QCD results are interpreted in fixed order without any attempt 
to improve the formulae through the inclusion of guesses for higher order coefficients. An estimate 
of the scale dependence is easily obtained through the evaluation of a variant of eq. (pi]) where 
i?NS is calculated for a general t'Hooft scale ff 2 . Adopting a s (10.5 GeV) = 0.177 corresponding to 
a s (Mz) = 0.12 and varying /j, 2 between s/4 and 4s the predicted value of R varies by -2 and +0.2 
permille. This is well below the anticipated experimental precision. Alternatively we may include 
in eq. an term with the coefficient based on a recent estimate in |l(|. This would lead to a 
decrease in R by 1.7 permille, again far below the forseeable experimental accuracy. 

Let us now discuss the impact of initial state radiation, the running a and the T resonances 
on the hadronic cross section, as observed at around 10 GeV under realistic experimental conditions. 
Lower energies contribute again through initial state radiation (eq. If we restrict the cutoff zq 

to a value of 0.25 which corresponds to a cut on the mass of the hadronic system of around 5 GeV, 
eq. ( |l~6|) for i?had can be applied also for o-q(zs') which appears in the integrand of (||). At the same 
time this cutoff excludes the region of the broad charmonium resonances which have not been well 
explored up today. A cutoff around z = 0.7 reduces the initial state radiation corrections further 
and eliminates contributions from the lower energy range completely. The vacuum polarization of 
the virtual photon has been discussed before for the case of the r lepton production and is identical 
for the hadronic cross section. The interference between the T resonances and the virtual photon 
is completely analogous to the leptonic case and leads to the same relative corrections. However, 



Table 4: Relative contributions from the T resonances to the hadronic cross section 
a(e + e~ — ► hadrons) for the energies E\ = 10.52 GeV and E 2 = 9.98 GeV and the two different 
values of the cutoff m m [ n = 5 GeV and 9 GeV. Interference terms and radiative tails are listed sep- 
arately. Also given are the continuum contributions and the resulting predictions for the total cross 
section. 





T(l) 


T(2) 


T(3) 


T(4) 


T(5) 


T(6) 


Ei = 10.52 GeV, m min = 5 GeV: 


Continuum: 


3.2228 nb 


Interf . : 


2.9 • 10~ 4 


3.1 • 10~ 4 


7.7- 10- 4 


-1.06- 10" 3 


-3.1 • 10" 4 


-9.7- 10" 5 


interf.: 


-1.0- 10" 4 


Rad. tails: 


5.88 • 10~ 3 


5.38 • 10~ 3 


1.212 • 10" 2 








r &d. tails: 


2.343 • 10~ 2 


Total: 


3.2980 nb 


E x = 10.52 GeV, m min = 9 GeV: 


Continuum: 


2.8515 nb 


Interf. : 


4.0 • 10" 4 


3.7 • 10" 4 


8.9- 10" 4 


-1.19- 10" 3 


-3.3 • 10" 4 


-1.1 • 10" 4 


interf.: 


3.9 • 10" 5 


Rad. tails: 


6.64 • 10~ 3 


6.08 • 10~ 3 


1.370 • 10" 2 


5 • 10" 5 






rad. tails: 


2.648 • 10~ 2 


Total: 


2.9272 nb 


E 2 = 9.98 GeV, m min = 5 GeV: 


Continuum: 


3.5568 nb 


Interf. : 


6.7- 10" 4 


-3.58 • 10" 3 


-4.6 • 10" 4 


-1.6- 10" 4 


-1.6- 10" 4 


-5.8 • 10" 5 


interf.: 


-3.74 • 10" 3 


Rad. tails: 


1.145 • 10" 2 


9 • 10" 5 










rad. tails: 


1.154- 10" 2 


Total: 


3.5845 nb 


E 2 = 9.98 GeV, m min = 9 GeV: 


Continuum: 


3.0666 nb 


Interf. : 


8.8 • 10" 4 


-4.1 • 10" 3 


-5.0 • 10" 4 


-1.7- 10" 4 


-1.6- 10" 4 


-6.1 • 10" 5 


J2 interf.: 


-4.14 • 10" 3 


Rad. tails: 


1.328 • 10" 2 


1.04 • 10" 4 










r &d. tails: 


1.338 • 10" 2 


Total: 


3.0949 nb 



the contributions from the radiative tails of the Breit-Wigner amplitudes are distinctively different. 
Without initial state radiation the tail from one resonance is given by 

_ i2^r e r had M 4 m 2 

°"had,R — ^ (s — M 2 ) 2 + T 2 M 2 ' 1 ' 

For the T(4, 5, 6) resonances an energy dependent width T(s) has to be used in (|l7|), which varies 
rapidly in the region of interest. For s below the BB threshold this leads to a drastic suppression of 



a(e + e" -> hadrons, m. , > m . ) [nb] 



had miir 



3.5 



2.5 



E 2 = 9.98 GeV 



E l = 10.52 GeV 



10 



in. 



[GeV] 



Figure 4: Dependence of the hadronic cross section a(e + e — > hadrons) on the cutoff in the minimal 
invariant mass of the hadronic system, m m i n; for the two energies E\ = 10.52 GeV, E 2 = 9.98 GeV 
and a s (M z ) = 0.120. 



these contributions. In the narrow width approximation eq. ( |i7|) implies 
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The contributions of interference terms and radiative tails from the individual resonances to the 
hadronic cross section are listed in Table |I] for the energies E\ = 10.52 GeV and E2 = 9.98 GeV and 



«s(Mf) 



0.120. Also given are the continuum contribution and the sum. The cutoffs m r 



5 GeV 



and 9 GeV are adopted. These correspond to the minimal and maximal values recommended for 
a precision measurement of i?had- The relative strength of the interference terms is identical for 
hadronic and leptonic final states. Again one observes the accidental cancellation at 10.52 GeV and 
the dominance of the (negative) T(2) contribution at the energy E 2 = 9.98 GeV just below this 
resonance. Compared to the leptonic cross section (Table |3|) the radiative tails are significantly more 
prominent — a consequence of the relatively large direct, non-QED-mediated hadronic decay rates of 
the T resonances. This leads to a positive correction of 2.3% and 1.2% for E\ and E 2 , if m m ; n = 5 GeV 
and to 2.7% and 1.3% if m min — 9 GeV, respectively. The difference in the strength of the T tails 
leads to an apparent variation of i?had around 1.3%. The cutoff dependence of the cross section is 
illustrated in Fig. ^. The strong dependence of the cross section for tight cuts is again clearly visible, 
suggesting a cut between about 5 and 9 GeV for E\ and 5 and 8.5 GeV for i?2- 



4. Uncertainties 



From the comparison of different radiator functions for initial state radiation the theoretical uncer- 



tainty from this source can be estimated to be below one permille. The error induced through the 
present simplified treatment of the hadronic vacuum polarization in the timelike region is estimated 
around two permille and could easily be reduced even further, if required. The combined theoreti- 
cal uncertainty from these and other effects is generously estimated below five permille. In addition 
photons from initial state radiation might be included in the invariant mass of the hadronic system 
or, conversely, photons from final state radiation may escape the detection. This 0(a) effect can 
only be evaluated for the concrete experimental analysis with the help of a Monte Carlo simulation. 
These theoretical uncertainties are significantly below the expected experimental error of roughly two 
percent, which is dominated by systematical uncertainties pj. 



5. Summary and Conclusions 

Precise predictions have been presented for the total cross section in the energy region explored 
presently by the CLEO experiment and at a future B meson factory. The present sample of nearly 
one million hadronic events allows for a small statistical error. These measurements will determine 
the value for a s under particularly clean conditions similar to the Z line shape measurements but at 
a different energy. When compared to the experimental results from Z decays, a determination of R 
with a precision of 2.5% would evidently demonstrate the running of a s between 90 and 10 GeV. A 
precision of 0.3% would be competitive with the a s measurements from the Z line shape which are 
based on the combined results of all four LEP experiments. 
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